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In this paper, we discuss a simple deterministic model for 
a field driven, thermostated random walk that is constructed 
by a suitable generalization of a multi-baker map. The map is 
a usual multi-baker, but perturbed by a thermostated exter- 
nal field that has many of the properties of the fields used in 
systems with Gaussian thermostats. For small values of the 
driving field, the map is hyperbolic and has a unique SRB 
measure that we determine analytically to first order in the 
field parameter. We then compute the positive and nega- 
tive Lyapunov exponents to second order and discuss their 
relation to the transport properties. For higher values of the 
parameter, this system becomes non-hyperbolic and posseses 
an attractive fixed point. 

05.45 + b, 05.70 + Ln 



I. INTRODUCTION 

In the past several years a great deal of attention has 
been devoted to computer and analytic studies of the 
chaotic properties of fluid systems subjected to external 
fields and to Gaussian thermostats which maintain a con- 
stant kinetic or total energy in the system, in the pres- 
ence of the field, Jl|,||. The interest in this subject stems 
not only from the method's value as a means of simu- 
lating non-equilibrium flows and computing their prop- 
erties, but also because there is a connection between 
transport properties, non-equilibrium fluctuations, and 
the underlying microscopically chaotic properties of the 
fluid. This connection has been explored from computa- 
tional ||-|| and analytic |]Hl|| points of view. The pur- 
pose of this paper is to describe a model system in which 
the transport and dynamics of a thermostated system can 
be studied in great detail, and in which one can explicitly 
construct the SRB measure [[id] and describe such prop- 
erties as the transition from hyperbolic to non- hyperbolic 
behavior, and related phenomena. These properties have 
been explored in previous work [ ^0|Jl5| , but have not yet 
been studied in great detail, due either to the complex- 
ity or to the simplicity of the models treated up till now 
p| ^3| . The model discussed here allows one to gain 
some insights into the general class of properties of ther- 
mostated systems, while keeping the analytical and com- 
putational difficulties to manageable proportions. It is 
one of the few cases known so far where one can check 



some of the general properties of thermostated systems 
on a specific model. 

The model we consider is a variant of the multi-baker 
maps studied by Gaspard and cowor kers [HUH, which 
are deterministic models for the diffusion of a particle on 
a one-dimensional lattice. The map considered here has, 
in addition, an external driving field which is constructed 
so as to model the effect of a thermostated electric field 
on charged particles in a two dimensional setting. We 
present the model and then calculate the chaotic prop- 
erties at small values of the external field. We obtain an 
expression for the stationary state SRB measure to first 
order in the applied field, and the positive and negative 
Lyapunov exponents to second order in the applied field. 
This allows us to verify the interesting relations between 
the rate of entropy production, the zero field diffusion 
coefficient, the drift velocity, and the sum of the Lya- 
punov exponents 

§13,011]. We conclude with a brief 
discussion of the transition to non-hyperbolic behavior 
as the field increases beyond a certain value, and discuss 
the connection of our model to other types of field-driven 
multi-baker maps JllHL 



II. THE NON-LINEAR MULTI-BAKER MAP 

We begin by considering a simple multibaker map that 
acts on the (x, y)-coordinates of particles, and that mod- 
els a random walk on a one-dimensional lattice of unit 
spacing. The map, defined on Z x [0, l] 2 , replaces the 
n,x,y coordinates of a particle by Mo(n, x, y) 



y 

M (n,x,y) = (n - 1,2a;, -), 



< x < 1/2, 



= (n + l,2i-l,^yi), 1/2<x<1. (1) 

Here n represents the position of the random walker on 
the line and x, y can be seen as bookkeeping variables 
keeping track of the deterministic cause of the apparent 
random walk. The subscript, zero, on M indicates that 
this is the map defined without an external field, which 
we introduce shortly. The map Mq is time reversal sym- 
metric. That is, there exists an involution operator, T, 
which acts on the x, y- variables, but not on the box in- 
dex n, and is given by T(x,y) = (1 — y, 1 — x), with 
T 2 = 1, where 1 is the identity operator in R 2 , and such 
that T o Mq o T(n,x,y) = M " 1 (n, x, y). If we consider 



1 



periodic boundary conditions, the invariant measure is 
uniform and has Lyapunov exponents 



A (o) = _ A (o) =ln2 _ 



(2) 



Next we suppose that the particles are also acted upon 
by a thermostated electric field, whose action we now 
model. Our final map will then be a composition (to be 
explained below) of the field map (|^) with the multi-baker 
map (0). 

The modelling of the field map can be done by con- 
sidering the action of a thermostated electric field on 
a (continuously) moving particle, where the thermostat 
maintains a constant kinetic energy for the particle. The 
equation of motion of a particle in such a field is given 

by [01 

dp $ (E-p\^ 
di =qE - q [—) P - 

If 9 is the angle that the particle makes with respect to 
the direction of the electric field, then 9 changes in time 
as 9 = —qEsm9/p, with solution 



,'(9(0) \ / qE 
9(t) = 2 arctan [ tan [ - I exp I 1 



A time discretized version of this equation is obtained by 
defining an angle at discrete times 9 n = Trx n , and letting 
x n satisfy 



2 

7T 



x n+ i = — arctan ^tan ^ " ) < 



Here x n £ [0, 1] and a = ^p-T, with r the time step, 
which for the time being we set equal to unity. Note that 
we restricted our attention to angles 9 E [0,ir], taking 
advantage of the symmetry 9 «-> —9. 

We now introduce a one parameter family of maps of 
the unit interval onto itself by 



<fia ( x ) = — arctan 

7T 



( ( 
I tan I ■ 



2 J 



(3) 



We note the following property : 

ip a (x) = 1 - <P-cc(l - X), (4) 

which implies the time reversal symmetry 

To (<p a ,l) oT(x,y) = (x,^- 1 ^)), (5) 

where T is defined above, and the operator (ip a , 1) acts 
on a point {x,y), say to produce (<p a (x), y), and we have 
denoted the identity operator acting on the y-coordinate 
by 1. We also point out the fact that the map (p a {x) has 
the property, under successive iterations, that 



which shows that the field map is a good discretization 
of a continuous process in time. 

We can now construct the following non-linear (or field 
driven) multi-baker map as a time reversal symmetric 
composition of the multi-baker map (0) with the field 
map (|) : 

M a (n, x, y) = (1, 1, ip a ) o M a ° (1, <p a , l)(n, x, y), (7) 
which takes the explicit form 

y 

M a (n,x,y) = [n - 1, 2<p a (x), <p a (-)), 

< X < <p- a {^), 

y + 1 

= (n + 1, 2tp a (x) - 1, <p a (— — ))) 

<P-«(\) < ^ < 1. (8) 

The leftmost identity operators in each of the field maps 
in (0) act on the cell index, n, and express the fact that 
the field maps do not change the value of the cell in- 
dex. Only the baker map moves points from one cell 
to the next. We refer to Fig. (5) at the end of this 
paper for an illustration of the projection of M a along 
the s-interval. The time reversal symmetry of this map 
T o M a o T(x, y) — M~ 1 {x 1 y), with T defined as above, 
follows straightforwardly from (^, ||) . 



III. HYPERBOLIC REGIME 

As long as a < ln(2), M a is expanding along the x- 
direction, i.e. 



dM a: 
dx 



> 1, 



(6) 



so that standard theorems guarantee the existence and 
uniqueness of an SRB measure fl9| . In this section we 
solve for the invariant density and give an analytic ex- 
pression of this invariant measure. 

We want to find the stationary eigenfunction, equiv- 
alently, the invariant density, p(n,x,y), of the Perron- 
Frobenius operator for a system with periodic boundary 
conditions. This implies that p does not depend upon n, 
but only on x,y and satisfies the equation 

x x 
P{x,y) = ip'_ a (~)v'_ a (y)p(<p- a (-),2ip- a (y)), 

0<y< ¥>a(i), 

M\)<y<^ (9) 

where the prime denotes the derivation with respect to 
the argument. We solve this equation by expanding the 
density in powers of the field parameter, 
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p(x,y) ~ 1 + ap {1 \x, y) + a 2 p (2) (a;, y) + o(a 3 ). (10) 

For small a, the low field regime, we expand tp a in 
powers of the field parameter : 

a a 2 
<p a (x) ~ x sin(7rx) + — sin(27ra;) + o(a 3 ). (11) 

The first order correction to the invariant density is found 
by decomposing it in Fourier modes : 



P (1) {x,y) =^2[a k (y)co$(2Trkx) + b k (y)sin(2Trkx)]. (12) 

fc=0 



Inserting the expansion for p, Eq. fll0|), and Eq. 
in Eq. (||), we can find the a^'s and b k s. However, as 
a result of the phase space contraction, the density is a 
singular function of the y coordinate so that we cannot 
represent the a k s and b k 's in terms of standard functions 
p0| . For our purposes, it is enough to perform a partial 
integration of the density along the y direction so as to 
obtain continuous coefficients for the Fourier modes. We 
thus define 

MV)= [ V a k (y')dy'; B k (y) = [ % ' b k (y')dy' '. (13) 
Jo Jo 

which, to lowest order in a, are found to satisfy the re- 
cursion relations 



Mv) 



-A {2y) + — + -sm{Tvy) 

I IT IT 



7T ^-^ 



k' odd 



B k >(2y) 
k' ' 



1 



A (2y-1) 

7T 



k' odd 



0<y<~, 
2(1 -v) 1 . / x 

+ — — + - sin(Try) 

7T 7T 

B k ,{\) B k ,(2y-1) 
k' 



k' 



- 2 <y<i, 



(14) 



A k {y) = \A 2k {2y)-- nT ^ 



?r(4r - 1/4) 



E 



k 1 



B k >(2y), 



0<y< 



1 



2 

(A, t {i) + A M {2y-\))- -j-^iL_ 
§ E S ^( B *'d)-B t '(2 ! /-l)), 

A;' odd 

^ < y < i, (15) 



1 4fcw 
B ^) = 2 B ^) + ff(4fca _ y i/4) 



4fc x 
+ — > 

/c' odd 



iF^^' (2j/) ' 



0<y < 



1 

2' 
4fcy 



+^(^-l)) + - 7r(4fe2 _ 1/4) 



4fc x ^ 1 

+ ~ ,,,, 4fc 2 - fc' 2 

fc' odd 



(A fe ,(l)-^(2y-l)), 
1 



<y < l. 



(16) 



In particular, for y — 1, we find A k (l) — A2 k (l) and 



B k (l) = B 2k (l) 
spcctively, 



4fc 



?r(4r - 1/4) 



whose solutions are, re- 



A k (l) = 0; S fe (l) = -]T 



2™fc 



7r^ 2 2 (" +1 )fc 2 -l/4 



(17) 



Therefore the projection of p^'(x,y) along the IE- 
direction, is a smooth function of x given by 

„1 oo 

p {l \x)= / p( 1 \x,y)dy = VB fe (l)sin(27rte). (18) 
•'o fe=l 

To compute the coefficients (y) and 73 & (y) numerically, 
we need to select a cut-off value k max of k beyond which 
we set all the coefficients to be zero. In Figs. (l)-(3), 
we show Aq and the first five A k 's and B k 's, respectively, 
computed by setting k max = 250. 

Although we have found the invariant density only to 
first order, we can now compute the corrections to the 
Lyapunov exponents to second order in a. For A+, we 
obtain 



i) 



i ,i 

dx / dyp(x,y)1n(2<p a (x)) 
lo 



= ln(2) - a 2 f ~ + J dxp {1) (x) cos(nx) 

kB k {\)\ 



M2)-W7 + rE 



A 



(o) 



a 2 (0.6661513± 10~ 7 ). 



(19) 



Here we used the normalization condition for p(x, y) 
which requires that 



l ,i 
dx / dyp^ (x, y) = 0. 
n Jo 



(20) 



To compute A_, we need the full expression of p^', 
Eqs. : 
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A_ = / dx I dyp{x,y)\n(-(p' a (y/2)) 

Jo Jo 1 

+ ' 



dx I dyp{x,y)\jx{\v' a {{y + l)/2)) 

(p- a (l/2) Jo L 



<P- a (l/2) -M2) 



2a a 2 
it 4 



1/2 



<r I dx I dypi(x,y)cos(Try/2) 
h Jo 



+ (l-^_ a (l/2)) -ln(2) + 



2a a 2 

TT 4 



+ a 2 dx dypi(x,y)sm(iry/2) 

Jl/2 Jo 

= A (0)_ a2 A + ± + I E Mll 

\ 4 7T 7T ^ fc 

\ odd 

+ ~y dyA (y){sin(iry/2) + cos(iry/2)) 
+ 1 /^E ^(sin(^/2)-cos(^/2))) 

fcodd K I 



(21) 



It is not straightforward to compute these integrals nu- 
merically because of the irregularity of the functions ( [l4| - 
|l6| ) and the number of different terms involved in their 
expressions. We can nevertheless estimate ( ^T|) within 
some good accuracy. To this purpose, we proceed by a 
number of algebraic manipulations. 

We first substitute for B k (y) the expression 



Bk(y) = yB k (l) + kf k (y), 



(22) 



where the functions fk(y) are found to satisfy the rela- 
tions 

My) = /a*(2y) + -Y, ,,2 1 7jS A v<. 2 v)> ( 23 ) 

k' odd 

for < y < 1/2, and fk(y) is odd with respect to 1/2, 
viz. f k {y) = -/fc(l - y). 
In terms of these, we have 



A k {y) = \A 2k {2y)-- ]T 



J2 



tt f-^4fc 2 -fc' 2 

k' odd 



/ fe '(2y) 



TT I 16fc 2 - 1 + E 4fc 2 - /c' 2 

\ k' odd 



4fc' 



5^(1) , (24) 



for < y < 1/2, and A k (y) is even with respect to 1/2, 
viz. A k (y) = A k (l - y). 

We can now find an upper bound on the magnitude of 
f k (y). Indeed, A k (y) is everywhere negative and is mini- 
mal at y = 1/2 (see Fig. (2)). Thus 



i^dii^+E i w^ Bk ' {1) )- (25) 

\ k' odd / 



Now, f k (y) is negative between and 1/2 and reaches its 
minimum at y — 1/4. Thus, from Eqs. ( p3| , p5| ), 



IM</, 



2 x-^ 1 

'max — X 

^ ~ ^ J-i A 4fc 2 - 

Ar odd 



4fc' \ 
2-772^(1) • (26) 



4fc 

Next, we rewrite Ao(y), Eq. (p^), in terms of its 
Fourier modes : 

My) = ~f]Gk cos(2irky) + H k sin(27rfcy). (27) 
We find 



fc=0 



A:' odd 



fc' odd 



gfc^l) 



?r(4fc 2 - 1) 
+ 2 E / d Vfk'{y)cos{kiTy), fcodd, 

fc' odd 1 ' 



1 



Gfc ~ 2 Gfe / 2 " ^(4fc 2 -l)- 
H k = 0, k odd, 



A: even, 



^ E ^ + 2 E t dyfAy)sHknyl 

fc' odd fe' odd 170 

A even, (28) 
Therefore and with the help of Eqs.(p5|, p6| ), 



k odd 

x y k odd 



-(sin(7ry/2) — cos(iry/2)) 
B k (l) 



[ d yfk(y)(sHny/2) ~ cos(7n//2)) 



/c odd 

2 /4 



A; odd 



max 
k 



(29) 



and 



dyA (y){sin(iry/2) + cos(?n//2)) 



4G 4 y> G k 

7T 2 7T 2 ^ 16/C 2 - 1 

k— 1 

4Gq ^ J_ 

TT 2 ^ 2 ^ \ ^ 2 n 2 2n 

k odd \n=0 



fc 2 -l 



4 



oo n — 1 



' i ' i 2^ 2 2n 



71 = 1 j=Q 



i+4 1,2 



^ * 2 , V n 2 " 2 2n+4 fc 2 - 1 

odd \n— 

oo n — 1 
\ ~> \ ^ 1 ff2"~Jfc 
+ oi 9 2n+4 .2 _ 



n=l j=0 



n(± yf 1 1 

+ V ^ A, ^ 2"" 1 2 2 ™ +4 fc 2 - 1 
\ odd n— 



max 



2 V"^ fk 

k' odd 



where we introduced the notations 



G 



and 



(0) 



21.2 



TT Z k 



ff2"-J/c 



fe' odd 



-By(l) 
k' 



7r(4fc 2 - 1) 



7r(2 2(n- J + D A .2_ 1) - 



(30) 



(31) 



(32) 



Grouping Eqs.(|29|, |30|) together with Eq.(|21|), we can 

give an estimate of the second order correction to the neg- 
ative lyapunov exponent by performing straightforward 
numerical summations : we find 



A_ = A 



(0) 



a 2 (1.9937± 0(0.04)). 



(33) 



In Fig. (4) we compare the values (19) and (|33j) of the 
second order corrections to the Lyapunov exponents to 
numerically computed ones. 

In the next section we take the macroscopic limit and 
relate the sum of the Lyapunov exponents to the drift 
velocity and the zero-field diffusion coefficient. 



IV. MACROSCOPIC LIMIT 

In order to take the macroscopic limit of this diffusive 
process, we consider the function Wt(n) which is defined 
to be the total probability of finding a particle in cell n 
at time t, and is given by 



W t {n) 



dx / dyp t (n,x,y), 
o Jo 



(34) 



where pt(n, x, y) is the time-dependent density of the ex- 
tended version of M a and whose time evolution is derived 
by generalizing Eq. (^) in a straightforward way. Once 
we have eliminated the internal variables, (x,y), we no 
longer have a deterministic process, but instead have a 
random process which should, in a macroscopic limit, be 
described by a suitable Fokker-Planck equation with a 



drift term, representing the effect of the external field. 
The macroscopic limit is taken by scaling the space and 
time parameters, and then taking an appropriate scaling 
limit. This procedure was described for models of this 
type by Tel, Vollmer, and Breymann [ [TTIJT^ ], and here 
we simply outline the process. 

We have previously introduced the time step r, and 
we replace the times t and t + 1 in the Frobenius-Perron 
equation by Tt and Tr + r, respectively, where T » 1. 
Similarly, we scale the length of the elementary cells, n, 
by making them have a length a on a side. We then 
replace n and n + 1 in the Frobenius-Perron equation by 
Na and Na + a, respectively, where N » 1. Also the x 
and y variables have to be scaled by the factor a, as well. 

We can now proceed to the derivation of the Fokker- 
Planck equation from the Frobenius-Perron equation for 
Pt(n,x,y). We obtain a Fokker-Planck equation for a 
field-driven random walk by considering the difference 
W TT+T (Na) - W TT (Na). 

Using Eqs. we find 

rv- a (i/2) 

Wtt+t(Nci) = / dxpT T (Na + a,x) 

Jo 

dxpTr(Na — a, x), 

where pT T (Na,x) — a dypTr(Na, x,y). We have 
scaled x and y so that their values are in the interval 
< x, y < 1. Therefore, 

Wtt+ANci) - W T (Na) = 
rv- Q (i/2) 

dx(pT T (Na + a, x) — pTr{Na, x)) 



dx(pT T (Na — a,x) — pTr(Na,x)). (35) 

V-«(l/2) 

Expanding pT T {Na ± a,x) about pT T {Na,x) and 
Wtt+t{No) about Wtt{Ncl) respectively, and introduc- 
ing the coordinate X — Na, and the time t = Tt, we 
get 



dWt(X) 




dxp t {X, x) 
dxpt(X, x) 1 + a 



2 1 d 2 W t (X) 



¥>-«(l/2) ) 2T dX 1 

The drift velocity is given by, 



v t {X) = ~x 

T 



(36) 
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d 

dX 



r<P- a (l/2) pi 

/ dxp t (X, x) + / dxp t (X,x) 



d f 1 

dxj tort*'*) 



(37) 



To proceed further, we need to take the macroscopic limit 

2 

where a — > 0, t — > 0, t^t = D and v is finite and nonzero. 
This implies that the electric field, E, becomes infinite 
as a~ 2 . 

For long times, t, we are going to replace the drift 
velocity in the Fokker-Planck equation by its stationary 
state value. Suppose t is large enough that we are nearing 
a steady state. Then in the limit of small a and r, we 
may write the solution of the Frobenius-Perron equation 
p t (X,x) in ( |37| ) as Wt(X)p(x) where p(x) is determined 
by an equation easily obtained by integating the equation 
@ over y (|Ts|), and W t (X) is close to, but not quite, 
a constant. Then the X dependence drops out in the 
expression for the drift velocity, and we find that the 
stationary state drift velocity is simply (assuming that 
the density is normalized to a unit cell) 



_ o 

v = — 

T 



¥>-«(l/2) 



dxp{x) 



<p_„(l/2) 



dxp{x) . (38) 



Eq. (|3^) is the Fokker-Planck equation corresponding 
to a stochastic diffusive system with a drift. In the case 
of periodic boundary conditions that we consider here, 
we know from thermodynamics that the rate of entropy 
production in the stationary state is due solely to the 
existence of a current driven by the external field and for 
which the rate of entropy production is given by 



v 

d' 



(39) 



where the zero-field diffusion coefficient for this process is 

2 

D = S— . We mention that in the limit of zero field, and 
periodic boundary conditions, the entropy production in 
the stationary state vanishes, since in this limit, there 
is no steady state drift, and the distribution function 
p{X, x) is constant both in X and in x. 

With the stationary state density (|18| ) computed in the 
previous section, the drift velocity (p8f) is found to be 



2aa 



ITT 

act 



\ k odd / 



(1.15217813). 



(40) 



Hence, the entropy production rate is, Eqs. (|39j- 140|) , 

a 2 

cr = 2.65502889—. (41) 

T 

According to the usual arguments for thermostated sys- 
tems, one expects that the rate of phase space contraction 



given by the negative of the sum of Lyapunov exponents 
should be equal to the macroscopic rate of entropy pro- 
duction [ |To| , p~3| ] . For our case, it is possible to verify this 
relation analytically, since we have been able to calculate 
all of the relevant quantities. The phase space contrac- 
tion rate is given by 



(A+ + A_) = (2.66 ±0(0.04))- 



(42) 



and Eqs. (|d],[l2]) give consistant values. Thus, our field 
driven random walk model has a well behaved macro- 
scopic limit, provides an example of the correspondence 
between the macroscopic and microscopic relations for 
entropy production, and is analytically tractable. 



V. CONCLUSIONS AND DISCUSSION 

In this paper, we have shown that it is possible to con- 
struct a non-linear version of the multibaker map (i.e. 
it shares the topology and time-reversal symmetry of the 
original multi-baker but is not piecewise linear) that sim- 
ulates the action of an external field on a diffusive pro- 
cess. The field curves the branches of the map and is 
responsible for the phase space contraction that induces 
a stationary state on an attractor (it fills the whole phase 
space but its information dimension |14| is fractional as 
a consequence of the difference between the second order 
corrections to the Lyapunov exponents). 

One of the motivations of this work was to provide 
an analytically tractable map which shows some of the 
properties of the periodic Lorentz gas where the parti- 
cle moves among the scatterers in a thermostated elec- 
tic field. The structure of our non-linear baker map is 
sufficiently simple that we were able to compute analyt- 
ically the stationary state SRB measure using a pertur- 
bation expansion in the field parameter. This allowed 
us to compute the positive and negative Lyapunov expo- 
nents, whose values we showed are consistent with that 
of the drift velocity. We also were able to compute the 
irreversible entropy production and showed that it is in- 
deed given, apart from some factors, by the sum of the 
Lyapunov exponents. 

We remark that at values of the field parameter larger 
than ln(2), the map loses its hyperbolicity. This is illus- 
trated in Fig. (5). In fact (0,0) becomes an attractive 
fixed point of the reduced map, which means that, on the 
lattice, all particles eventually move ballistically around 
the ring towards decreasing n's. The case a = ln(2) is of 
particular interest. Indeed, the origin is an intermitent 
fixed point and, as a consequence, points can spend arbi- 
trary long times in its vicinity This can be seen to give 
rise to anomalous diffusion pi| , p2| . 

We regard the model given here as the simplest of a 
class of similar models which can be generated by vary- 
ing an additional parameter modelling a magnetic field. 
Elsewhere [E3i we will describe this class of models in 
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much more detail, because they show a wide variety of 
features both in the hyperbolic and non-hyperbolic re- 
gions, including sequences of period adding bifurcations. 

In a subsequent paper pif we will consider boundary 
conditions other than periodic. It is in particular an im- 
portant question to check whether, with flux boundary 
conditions, the chain sustains a stationary state with an 
almost linear gradient of density. As we showed in an- 
other paper jl3| , some other models of field driven multi- 
baker maps fail to have this behavior in a large system 
limit and, as a consequence, do not have physically rel- 
evant thermodynamics (the entropy production rate de- 
pends on the choice of the partition). 

Another interesting perspective will be to calculate non 
linear corrections to the diffusion coefficient so we can go 
beyond the linear response theory. As shown by other au- 
thors , relevant tools for this study are the zeta func- 
tions and Pollicott-Ruelle resonances. Computing the 
first order corrections to the eigenvalue spectrum of the 
Perron-Frobenius operator will be an important step to- 
wards understanding non-linear diffusion for this model. 
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FIG. 2. The first five Ak(y) computed with a cut-off value 

^max — 250. 



FIG. 4. A comparison of the values of the second order 
corrections of the Lyapuonv exponents computed from Eqs. 
([l9|, |il|) (solid lines) and by following a trajectory of 10000 
steps for 10 values of a ranging from to .1 (x's correspond 
to the positive Lyapunov exponent, +'s to the negative one). 
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